Single-flavor CSL phase in compact stars 
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Abstract. We suggest a scenario where the three light quark flavors are sequentially 
deconfined under increasing pressure in cold asymmetric nuclear matter as, e.g., in neutron 
stars. The basis for our analysis is a chiral quark matter model of Nambu-Jona-Lasinio 
(NJL) type with diquark pairing in the spin-1 single flavor (CSL), spin-0 two flavor (2SC) 
and three flavor (CFL) channels. We find that nucleon dissociation sets in at about the 
saturation density, no, when the down-quark Fermi sea is populated (d-quark dripline) 
due to the flavor asymmetry induced by /3-equilibrium and charge neutrality. At about 
3no u-quarks appear and a two-flavor color superconducting (2SC) phase is formed. The 
s-quark Fermi sea is populated only at still higher baryon density, when the quark chemical 
potential is of the order of the dynamically generated strange quark mass. We construct two 
different hybrid equations of state (EoS) using the Dirac-Brueckncr Hartree-Fock (DBHF) 
approach and the EoS by Shen et al. in the nuclear matter sector. The corresponding 
hybrid star sequences have maximum masses of, respectively, 2.1 and 2.0 M Q . Two- and 
three-flavor quark-matter phases exist only in gravitationally unstable hybrid star solutions 
in the DBHF case, while the Shen-based EoS produce stable configurations with a 2SC 
phase component in the core of massive stars. Nucleon dissociation due to d-quark drip at 
the crust-core boundary fulfills basic criteria for a deep crustal heating process which is 
required to explain superbusts as well as cooling of X-ray transients. 
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1. INTRODUCTION 



The phenomenology of compact stars is intimately connected to the EoS of matter 
at densities well beyond the nuclear saturation density, n = 0.16 fm -3 . Compact 
stars are therefore natural laboratories for the exploration of baryonic matter under 
extreme conditions, complementary to those created in terrestrial experiments with 
atomic nuclei and heavy-ion collisions. Recent results derived from observations of 
compact stars provide serious constraints on the nuclear EoS, see pj and references 
therein. A stiff EoS at high density is needed to explain the high compact-star 
masses, M ~ 2.0 M , reported for some low-mass X-ray binaries (LMXBs), e.g., 
4U 1636-536 [2|, and the large radius, R > 12 km, of the isolated neutron star 
RX J1856.5-3754 (shorthand: RX J1856) Q. Another example is EXO 0748-676, 



an LMXB for which the compact-star mass and radius have been constrained to 
M > 2.10 ±0.28 M and R > 13.8 ±0.18 km [j. However, the status of the results 
for the latter object is unclear, because the gravitational redshift z = 0.35 observed 
in the X-ray burst spectra j^] has not been confirmed, despite numerous attempts. 
While compact-star phenomenology apparently points towards a stiff EoS at high 
density, heavy-ion collision data for kaon production p] and elliptic flow [3] set an 
upper limit on the stiffness of the EoS |l|. 

A key question regarding the structure of matter at high density is whether a 
phase transition to quark matter occurs inside compact stars, and whether it is 
accompanied by unambiguous observable signatures. It has been argued that the 
observation of a compact star with high mass and large radius, likewise reported 
for EXO 0748-676, would be incompatible with a quark core [4j, because the 
softening of the EoS due to quark deconfinement lowers the maximum mass and 
the radius in comparison to the nuclear matter case. However, Afford et al. [8] have 
demonstrated with a few counter examples that quark matter cannot be excluded 
by this argument. In particular, for a recently developed hybrid star EoS based 
on the DBHF approach in the nuclear sector and a three-flavor chiral quark model 
(lo| . stable hybrid stars in the mass range from 1.2 M Q up to 2.1 M q h ave been 
obtained, in accordance with modern mass-radius constraints, see also jljj] . In that 
model, a sufficiently low critical density for quark deconfinement was achieved 
with a strong diquark coupling, while a sufficient stiffness for a high maximum 
mass of the compact star sequence was obtained with a repulsive vector meanfield 
in the quark matter sector. The corresponding hybrid EoS for symmetric matter 
was shown to fulfill the constraints from elliptic flow data in heavy-ion collisions. 
In the present work we want to discuss a new scenario of quark deconfinement, 
which could play an important role in asymmetric matter, in particular for the 
phenomenology of compact stars. 

Chiral quark models of the NJL type with dynamical chiral symmetry breaking 
have the property that the restoration of this symmetry (and the related quark 
deconfinement) at zero temperature is flavor specific. When solving the gap and 
charge-neutrality equations selfconsistently one finds that the chiral symmetry 
restoration for a given flavor occurs when the chemical potential of that flavor 
reaches a critical value that is approximately equal to the dynamically generated 
quark mass, /i/ = /i c ~ rrif, where / = u,d,s. In asymmetric matter the quark 
chemical potentials are different. Consequently, the NJL model behavior suggests 
that the critical density of deconfinement is flavor dependent, see Fig. [IJ In this 
approach the down quark flavor is the first to drip out of nucleons as the density 
increases, followed by the up quark flavor and eventually also by strange quarks. 
This behavior is absent in simple and commonly applied thermodynamic bag model 
equations of state since they are essentially flavor blind. 

Under the /3-equilibrium condition in compact stars the chemical potentials of 
quarks and electrons are related by /x^ = fi s and ^ = /i u ±/i e . The mass difference 
between the strange and the light quark flavors m s 3> m u , rrid has two consequences: 
(1) the down and strange quark densities are different, so charge neutrality requires 
a finite electron density and, consequently, (2) \i d > fi u . When increasing the 
baryochemical potential, the d-quark chemical potential is therefore the first to 




FIGURE 1. Left panel: Chemical potentials of up and down quarks (strange quark sector not 
shown). With increasing quark chemical potential /j, q — + ^ti)/2 in isospin asymmetric matter 
the quark flavors pass sequentially the threshold (/i c ) for chiral symmetry restoration (deconfmc- 
ment), which entails nucleon dissociation. Right panel: Solution of the NJL gap equations under 
isospin asymmetry. 



reach the critical value /z c where the chiral symmetry gets (approximately) restored 
in a first-order transition and a finite density of d-quarks appears. Due to the finite 
value of n e , the u-quark chemical potential is still below fi c while the s-quark density 
is zero due to the high s-quark mass. A single-flavor d-quark phase therefore forms 
in co-existence with the positively charged nuclear-matter medium. 

Why has this interesting scenario been left unnoticed? On the one hand, bag 
models, which are commonly used to describe quark matter in compact star 
interiors cannot address sequential deconfinement. On the other hand, the single- 
flavor d-quark phase is negatively charged and cannot be neutralized in a purely 
leptonic background. This was a reason to disregard it in dynamical approaches 
like the NJL models. In the following we discuss the single-flavor phase for the first 
time under the natural assumption that the neutralizing background is nuclear 
matter. Since nucleons are bound states of quarks, the physical context in which 
such a mixed phase of nucleons and free d-quarks occurs is that of the dissociation 
of nucleonic bound states of quarks (Mott effect). 



2. PHASE TRANSITION TO QUARK MATTER: 
NUCLEON DISSOCIATION 



The task to develop a unified description of the phase transition from nuclear 
matter to quark matter on the quark level, as a dissociation of three-quark bound 
states into their constituents in the spirit of a Mott transition has not been 
solved yet. Only some aspects of such a description have been revealed within 
a nonrelativistic potential model 12, 13[ and within the NJL model 14]. We will 
consider a chemical equilibrium of the type n + n ^p + 3d, which results in a mixed 
phase of nucleons and down quarks once the d-quark chemical potential exceeds 
the critical value. This scenario is analogous to the dissociation of nuclear clusters 
in the crust of neutron stars (neutron dripline) and the effect may therefore be 
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FIGURE 2. Phase diagram in the plane of baryon and charge chemical potential. The dash- 
dotted line denote the border between oppositely charged phases. The nuclear matter EoS is 
DBHF (left panel) and Shen et al. (right panel). 



called the d-quark dripline. We approximate the quark and nucleoli components 
as subphases described by separate models. For t he nuclear matter subphase 

13, EE 01 with the 



we use two alternatives: (1) the DBHF approach [15|, [if 
relativistic Bonn A potential, where the nucleon selfenergies are based on a T- 
matrix obtained from the Bethe-Salpeter equation in the ladder approximation and 
(2) the EoS by Shen et al. j20|, which is based on a relativistic mean field theory 
and includes the contribution of heavy nuclei, described within the Thomas- Fermi 
approximation. The quark matter phase is described within a three-flavor NJL- 



21 . 122I . 23j |. The path-integral 



type model including diquark pairing channels [1 
representation of the partition function is given by 



T>qDq exp 




q(iqf — m + /ry ) q + C y 



Ant = Gs\Yl [(<? r a<?) 2 + (qil5T a q) 2 } + Vdo ^ 3do,a3do,a + Vdi 3dx3di \, ( 2 ) 

^ a=0 A=2,5,7 > 

where fi = + diagj(|, — |, — \)[iq + + Ag/is is the diagonal quark chemical 
potential matrix and m = dia.gf(m u ,m l j i ,m s ) is the current-quark mass matrix. 

For a = 0, To = (2/3) 1 / 2 !/, otherwise r a and A a are Gell-Mann matrices acting in 
flavor and color spaces, respectively. C = i r y 2r y is the charge conjugation operator 
and q = q^ . The scalar quark- antiquark current- current interaction is given 
explicitely and has coupling strength G$- The 3-momentum cutoff, A, is fixed 
by low-energy QCD phenomenology (see table I of (24J). The spin-0 and spin-1 
diquark currents are Jdo,a = q T iC^T A \Aq and j m = q T iC(^ji\ 7 + 72A5 + 73A 2 )g. 
While the relative coupling strengths 77^0 an d tJdi are essentially free parameters, 
we restrict the discussion to the Fierz values, 7]d = 3/4 and r\r>\ = 3/8, see (2o| . 
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FIGURE 3. Left panel: Compact star sequences. The phase structure of the core changes with 
increasing density, as indicated in the figure. Constraints on the mass come from 4U 1636 [2| and 
on the mass-radius relation from RX J1856 Q. Right panel: Hybrid equations of state used in 
the calculation of the compact star sequences. 



Color superconducting phases in QCD with one flavor have first been discussed 
in Refs. (26l . l27l . 28| . where also the special role of the spin-1 color-spin locking 
phase has been pointed out that here all quarks participate in the pairing with a 
gap of the order of 1 MeV or even below. This feature of the CSL phase is robust, 
as was demonstrated for our above isotropic ansatz for the spin-1 diquark current 
introduced in [2j|, and for its generalizations to the nonlocal case [30j and to a 
selfconsistent Dyson- Schwinger approach [311 ]. 

The gaps and the renormalized masses are determined by minimization of the 
mean-field thermodynamic potential, under the constraints of charge neutrality 



and /^-equilibrium. For further details, see [10|, |2l|, l22|, [23J]. In Fig. [2] we plot 
the thermodynamically favored phase in the plane of baryon and charge chemical 
potentials. 

The hybrid EoS corresponds to the dash-dotted lines in Fig. [2] and are con- 
structed such that the mixture of nuclear matter, quark matter and leptons is 
charge neutral. Using the hybrid EoS we calculate the corresponding compact star 
sequences by solving the Oppenheimer-Volkoff equations for hydrostatic equilib- 
rium. The hybrid-star sequences fulfill all modern constraints on the mass-radius 
relationship, see Fig. [3j For the DBHF hybrid EoS all stars with DBHF+CSL mat- 
ter in the core are stable equilibrium solutions, while the appearance of u-quarks 
and the associated formation of a 2SC subphase renders the sequence unstable. 
The situation is somewhat different for the Shen hybrid EoS, because in addition 
to Shen+CSL stars there are stable solutions with 2SC+CSL matter in the core. 
In both cases configurations with strange quarks in the core are unstable. 

The hybrid star sequences 'masquerade' as neutron stars [32j, because the 
mechanical properties are similar to those of nuclear matter stars and the transition 
from nuclear matter to the mixed phase is associated with a relatively small 



discontinuity in the density. Unmasking neutron star interiors might therefore 
require observables based on transport properties, which could be strongly modified 
in presence of color superconductivity. It has been suggested to base such tests of 
the structure of matter at high density on an analysis of the cooling behavior 



33l . l34j . l35l . |36| or the stability of rapidly spinning stars against r-modes [37|, [38 



It has turned out that these phenomena are sensitive to the details of color 
superconductivity in quark matter. 



3. BULK VISCOSITY AND URCA EMISSIVITY OF 
THE SINGLE-FLAVOR CSL PHASE 

According to [39] rotating compact stars would be unstable against r-modes in the 
absence of viscosity [40]. Constraints on the composition of compact-star interiors 



can therefore be obtained from observations of millisecond pulsars [371 . l38| . In such 
investigations the bulk viscosity is a key quantity. We therefore consider some 
relevant aspects here, starting with the two-flavor color superconducting phases 
following the approach described in Ref. [411 ] . 

The temperature-dependent bulk viscosity for the 2SC-CSL phase has been cal- 
culated self-consistently in j42| and is based on the flavor-changing weak processes 
of electron capture and (3 decay 

u + e~ — > d-\-u e , d — > u + e~ + v e . (3) 

It has been shown that the bulk viscosity is related to the direct URCA emissivity, 



which for normal quark matter was first calculated by Iwamoto [43| and can be 
expressed as 

914 7T 

£ ° ~ T680" G ^ VeVuVd T e 6 2 ue . (4) 

Here Gp is the weak coupling constant and ue is the angle between the up- 
quark and electron momenta, which is obtained from momentum conservation in 
the matrix element, see Fig. HI The triangle of momentum conservation holds for 
the late cooling stage, when the temperature is below 1 MeV and neutrinos are 
untrapped. Trigonometric relations are used to find an analytical expression for 
momentum conservation. To lowest order in 0d e the result is 

PF,d~PF,u ~PF,e - ~\vF,e 9%. ■ ( 5 ) 

For small angles 0d e — ue , so it is possible to obtain an expression for the matrix 
element of the direct URCA process. Following Iwamoto [43] one has to account 



either for quark-quark interactions to lowest order in the strong coupling constant, 
ots, & or the effect of finite masses (j7|): 



(6) 
(7) 



FIGURE 4. Direct Urea process in quark matter (left) and triangle of momentum 
conservation for it (right). 



From ©-(Cj) and the /3-equilibrium condition, /z^ = \i u + /i e , the angle 9^ that 
determines the emissivity (jl]) is obtained 




If interactions and masses are neglected, or the Fermi sea of one species is closed 
as in the single-flavor CSL phase, it follows that the triangle of momentum con- 
cervation in Fig. H|) degenerates to a line or can even not be closed. In that case 
the matrix element vanishes with the consequence that the direct URCA pro- 
cess does not occur, and also the bulk viscosity is zero. However, in the mixed 
nuclear-CSL phase there could be important friction and pair-breaking/formation 
processes, which we have not yet studied in detail. This could be an interesting 
issue for further investigation due to the large difference in the masses of baryons 
and deconfined quarks. 



4. MECHANISM FOR DEEP CRUSTAL HEATING 

Superbursts are rare, puzzling phenomena observed as a extremely long (4-14 
hours) and energetic (~ 10 42 erg) type-I X-ray bursts from LMXBs. They take 
place if the accreted hydrogen and helium at the surface burns in an unstable 
manner, which is the normal case [44]. As suggested by the authors of Ref. |4o| . 
superbursts could originate from accreting strange stars with a tiny crust and a 
core of three-flavor quark matter in the superconducting color flavor locked (CFL) 
phase, since it fulfills the constraints on matter properties from their superburst 
scenario. Of particular importance for the scenario is that in this phase the neutrino 
emissivity_and heat conductivity are suppressed by pairing gaps affecting all quark 



species [46l . 1471 . |48|. The mechanism underlying the superburst phenomenon is 
unstable thermonuclear burning of carbon in the crust at column depths of about 
10 e: cm" The carbon itself is a remnant of the burning of accreted hydrogen 



and helium at the surface. The ignition of observed superburst light curves takes 



place at a depth where the crust reaches temperatures around 6 x 10 8 K and column 
depths of about 10 12 g cm -2 . Such high temperatures in the crust at a certain depth 



are caused by deep crustal heating [50j, H3, The important ingredients for the 
strange star model of [45] are a thin baryonic crust of 100 to 400 m thickness, a 
sufficient energy release of 1 to 100 MeV per accreted nucleon by conversion into 
strange matter, a suppression of the fast direct URCA neutrino emissivity to the 
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and a thermal conductivity, k, of quark matter in the 
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FIGURE 5. Left panel: Density profiles of two stars with masses 1.4 M and 2.0 Mq. Note 
that the mixed phase of d-CSL quark matter with nuclear matter extends up to the crust-core 
boundary. Right panel: Energy release per nucleon as a function of the compact star mass. An 
upper estimate for the energy release from the conversion of DBHF nuclear matter to DBHF-CSL 
hybrid matter gives 12 (8) MeV for a compact star core with mass 1.4 (2.0) Mq. 

One of the main arguments for strange matter is the fact that superconducting 
phases, like the CFL phase, can suppress fast neutrino emission processes of all 
quark flavors and are able to fulfill the fusion ignition condition. However, as we 
have shown above, in the single flavor CSL-phase the fast direct URCA process is 
not possible at all, whereas slow neutrino cooling processes like bremsstrahlung of 
electrons and d-quarks occur. 

As one can see from Fig. [5] the energy release from the conversion of DBHF 
nuclear matter to DBHF-CSL hybrid matter gives 12 (8) MeV for a compact 
star core with mass 1.4 (2.0) M Q , which is in the range 1 - 100 MeV and could, 
in principle, explain burst ignition at appropriate depths for a suitable value of 
K. Therefore, the compact star does not necessarily need to be made of strange 
matter, but could be a hybrid star with quark matter in the d-CSL phase and 
a thin crust. Stejner et al. |44[ show that deep crustal heating mechanisms at 
the crust-core boundary like the conversion of baryonic matter to strange quark 
matter, which can fulfill the constraints of the superburst scenario do provide a 
consistent explanation of the cooling of soft X-ray transients too. Along the lines 
of this argument we claim that the d-quark drip effect at the crust core boundary, 
which leads to a mixture of nuclear matter with single-flavor quark matter in the 



CSL phase can serve as a deep crustal heating mechanism [53|, 54]]. Superbursts 
and the cooling of X-ray transients are not only consistent with quark matter in 
compact stars but may qualify as a signature! 



5. CONCLUSIONS 

In this contribution we have suggested a new quark-nuclear hybrid EoS for compact 
star applications that fulfills modern observational constraints from compact stars. 
Due to isospin asymmetry, down-quarks may "drip out" from nucleons and form 
a single-flavor color superconducting (CSL) phase that is mixed with nuclear 
matter already at the crust-core boundary in compact stars. The CSL phase has 
interesting cooling and transport properties that are in accordance with constraints 
from the thermal and rotational evolution of compact stars. It remains to be 
investigated whether this new compact star composition could lead to unambiguous 
observational consequences. We conjecture that the d-quark drip may serve as 
an effective deep crustal heating mechanism for the explanation of the puzzling 
superburst phenomenon and the cooling of X-ray transients. 
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